Reciprocal relationships in collective flights of homing pigeons 
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Collective motion of bird flocks can be explained via the hypothesis of many wrongs, and/or, 
a structured leadership mechanism. In pigeons, previous studies have shown that there is a well- 
defined hierarchical structure and certain specific individuals occupy more dominant positions — 
suggesting that leadership by the few individuals drives the behavior of the collective. Conversely, 
by analyzing the same data-sets, we uncover a more egalitarian mechanism. We show that both 
reciprocal relationships and a stratified hierarchical leadership are important and necessary in the 
collective movements of pigeon flocks. Rather than birds adopting either exclusive averaging or 
leadership strategies, our experimental results show that it is an integrated combination of both 
compromise and leadership which drives the group's movement decisions. 

PACS numbers: 89.75.Fb, 89.75.Kd, 89.75.Da 



I. INTRODUCTION 

The fast coherent movement of flocking birds is a 
fascinating phenomenon exhibiting apparent intelligence 
and coordination [1]. New monitoring technologies have 
meant that this collective behavior has recently attracted 
renewed interest from scientific and engineering commu- 
nities, allowing more conclusive analysis to be performed 
[2, 3]. It is a valuable topic whether all members of a 
flock contribute equally to the collective decision making 
and follow equivalent rules, or certain individuals have a 
greater influence on the decisions of the group [3-9] . 

The proposition that all contribute equally is some- 
times referred to as the "many wrongs" principle, and 
purports that individuals average their preferred direc- 
tions depending on interaction with their neighbors, lead- 
ing to a compromise in route choice [10, 11]. Conversely, 
the leadership hypothesis posits that one or a small num- 
ber of leaders are able to exert a disproportionate influ- 
ence on the group's movement decisions [12, 13]. Both 
theoretical and experimental arguments predict that the 
compromise of all members will make more accurate de- 
cisions than the leading of one or a small number in- 
dividuals, unless leaders have very different and superior 
information [4]. For example, homing performance of the 
pre-trained pigeons flying as a flock is significantly better 
than that of these birds released individually [14]. Yet, 
recent research has shown that time-varying hierarchi- 
cal decision making mechanisms do exist during pigeon 
flights [3, 6] — giving strength to the leadership hypoth- 
esis. 

Although new studies show that certain individuals in 
pigeon flocks are able to exert relatively more influence 
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on the movement decisions of the whole group, only the 
directed relationship (pointing from the leader to the fol- 
lower) has been studied [6, 7]. In addition to such di- 
rected links (representing a leader- follower relationship) 
selective coordinated behavior may also exist in collective 
motion, hence mutual relationships must also be exam- 
ined. For example, in the case of a perturbation caused 
by the terrain or a predator, the better strategy for birds 
would probably be to share the information of all group 
members to move rapidly to safety regardless of their in- 
dividual positions within a leadership hierarchy. Such 
mutual links represent a reciprocal relationship between 
a pair of pigeons, which appears nonrandomly in real-life 
directed networks [15, 16] and plays a significant role in 
the evolution of many biological systems [17, 18]. We em- 
phasize that the coexistence of compromise and leader- 
ship is not a contradiction, but a meaningful supplement 
to the hierarchical structure of pigeon fiocks. The main 
difference between the work in [6] and ours is as follows: 
they are concentrating on the leadership aspect; while we 
are focusing on local interactions (including both directed 
and reciprocal relationships). 

In this paper, we re-analyze experimental high- 
precision datasets of pigeon flocks to arrive at a more 
nuanced conclusion about the interactions and decisions 
in the collective dynamics of birds. Using quantitative 
methods from statistical physics [6], we find that both 
outcomes (directed and mutual links) coexist in the same 
flock flights. The mutual links represent a reciprocal re- 
lationship between individuals, which is a useful supple- 
ment to the well-defined hierarchical structure. Integrat- 
ing both directed and reciprocal links we uncover the 
complete topology of the network induced by the col- 
lective motion of a pigeon flock. Most significantly, our 
results imply that there is an integrated mechanism of 
decision-making in pigeon flocks: neither a leadership 
nor a compromise mechanism is clearly dominant, rather 
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both mechanisms coexist. 



perfect hierarchical structure induced by directed hnks 
in [6] may not be adequate to completely explain the 
collective behavior of pigeon flocks. 



II. RESULTS 



A. Reciprocal relationships by calculating pairwise 
correlations 



In the past tens of years, it is a very difficult mis- 
sion to explore the influence of individual members on 
a fast collective motion at all times. Recently, the ad- 
vance of GPS devices allow us to use sophisticated evalua- 
tion techniques to mine real flocking data [3] . Employing 
high-precision GPS in tracking pairs of pigeons. Biro et 
al. find that two birds compromise if they have less diver- 
sity on directional preferences, while either the pair split 
or one of them becomes the leader for a severe confiict 
[7]. Using lightweight GPS devices and analyzing data 
concerning leading roles in pairwise correlations, Nagy et 
al. show a well-defined hierarchy among pigeons belong- 
ing to the same flock [6]. In this study, we use the same 
datasets: 11 free flights and 4 homing flights. More de- 
tailed information on the datasets can be obtained from 
the website: http://hal.elte.hu/pigeonfiocks/. 

To investigate the influence that a given bird's behav- 
ior has on the other flock members, the temporal rela- 
tionship between the flight directions has been evaluated 
[3, 6]. The directional correlation for a pair of pigeons 

is Cij(r) — (^i(t) ■ Vj{t + t)'^ , where (• ■ • ) denotes time 

average and \^i{t) is the normalized velocity of bird i. 
When Cij{T) obtains its maximum value at the time de- 
lay Tij, Tij is called the optimal directional delay time. 
Negative r^- values mean that the ith bird falls behind 
the jth bird, which can thus be interpreted as a case of 
J leading. For each pair, we extract the positive value 
Tij = —Tji as a directed edge pointing from the leader to 
the follower. 

If individuals fly together in a flock, they will show 
a very similar velocity and a high correlation [19]. In a 
previous study, the authors study the leader-follower rela- 
tionship among pairs of pigeons whose directional corre- 
lation time delay is non-zero, and such links are directed 
[6] . The directed link indicates that a following bird tends 
to consistently copy the directional behavior of particular 
leading individuals. However, the directional correlation 
time delay may be near zero, which means a pair of pi- 
geons have a coordinated interaction with one another 
and there is a mutual (reciprocal) link between them. 
The frequency distribution of the directional correlation 
time delay for pairwise pigeons in all collective motions 
is shown in Fig. 1. Here we select 0.2 s as the time in- 
terval of time delay, for the sampling time interval of the 
original dataset and the resolution of time delay in the 
previous study [6] both are 0.2 s. The frequency of r = 
is the most frequent, meaning that there are typically 
many mutual relationships between the pairwise birds in 
pigeon flocks. Moreover, the result also implies that the 




FIG. 1: Histogram for the different values of the directional 
correlation delay time Tij between each pair of pigeons in 
all the flights. When the correlation coefficient reaches the 
maximum value Cij(r) at the time delay Tij, Tij is called 
the directional correlation delay time and Cij{T) is called the 
maximum correlation coefficient. Here we only consider the 
conditions of r,, > 0. 



Here we hope to explain why a correlation that de- 
cays from r = represents a "mutual link" instead of 
no interaction in our study. Actually, in some cases of 
collective motions, a decaying temporal correlation that 
is maximal at zero lag is an evidence that there is no in- 
teraction. For example, Katz and collaborators showed 
that the orientation correlation whose peak at zero time 
delay is significantly lower than the orientation correla- 
tion whose peak after zero [Fig. S8A in [20]]. On the 
contrary, in our study Cij (r) at t = tends to be higher 
than the values at the large time delay [Fig. 2], which 
means the strength of mutual links is stronger than that 
of directed links. 

Furthermore, we use multiple effective methods 
(Cij(T), motif, synchronization and multiscale analysis) 
to prove the existence of reciprocal interaction as follows. 
The values of the maximum correlation coefficient Cij (t) 
at the directional correlation time delay Tij have been 
shown in Fig. 2. The birds' movements under two con- 
ditions, free fiights and homing flights, are recorded [6]. 
The pigeon flock makes a circle-like route in free flights; 
while the group makes less direction turning during hom- 
ing flights [Fig. S6 in [6]]. Because of the more central- 
ized distribution of turning directions in homing flights, it 
is easier to obtain a larger Cij (r) , which results in the val- 
ues of Cij{T) at the same in homing flights [Fig. 2(b)] 
being larger than those in free flights [Fig. 2(a)]. In both 
cases, we can find that Cij(r) at the small Ty tends to 
be higher than the values at the large time delay. More- 
over, when the time delay is zero, the value of Cij(r) is 
the highest, which implies that these mutual links are 
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FIG. 2: The values of dj (r) decrease with Tij (s) for all the (a) 
free flights and (b) homing flights. The solid circles (squares) 
show the mean of Cijij) while the bars show the standard 
deviation of Cij{T). 



the most important relationship in the flock and more 
attention should be paid to such types of links. 

Although interaction and correlation are different and 
their relationship has been extensively discussed [19], us- 
ing the method of calculating pairwisc correlation [6], 
it is still difficult to determine whether the relationship 
between a pair of pigeons is a direct interaction or an 
indirect correlation. Here we develop a simple method to 
detect whether three links in a small subgraph (motif) 
are independent. A motif, defined as a small connected 
subgraph that recurs in a graph, is the basic functional 
unit of complex networks [21, 22]. A small motif with 
three nodes in the network induced from each flight are 
shown in Table I. 

TABLE I: Comparison of the values of the maximum corre- 
lation coefficient for the three links in the subgraph (motif) 
structure. Ci , C2 and C3 represent the maximum correlation 
coefficients for the three links, respectively, n is the whole 
number of the subgraph in all the flights. 

Motif n Ci > C2 Ci > C3 C2 > C3 ' ^ 

C2 > Ci 



173 54% 69% 64% 50% 



If the relationships among the three individuals are 
independent, the probabilities of Ci > C2, Ci > C3, 
and C2 > C3 should be about 50% respectively, while 
Ci > C3 along with C2 > C3 simultaneously should have 
a probability of 33%. However, we observe that all these 
relationships (except Ci > C2) occur significantly more 
frequent than we would expect [Table I], which means 
that B and C tend to have a lower correlation than that 
between A and B, and between A and C. Hence, the 



correlations among the three individuals are not consis- 
tent with three independent events and the directed link 
between B and C (C3) can be attributed to the other 
two high mutual correlations (Ci and C2). Again, our 
result implies that we can not neglect mutual links, or 
believe that the nodes with such type of relationship are 
no interaction. 



B. Multi-scale analysis of reciprocal links 

Although the pioneering researchers have obtained the 
detailed spatiotemporal data on the positions of indi- 
viduals during group movements, they do not use the 
data to study the dynamical variation of the relation- 
ship in collective motions [6]. To make clear the spa- 
tiotemporal relationship for a pair of mutual-link birds, 
we show the synchronization of B and G in the eleventh 
free flight [Fig. 3]. The high synchronization of two pi- 
geons again provides evidence that those pairs with mu- 
tual links {xij = 0) are strongly interacting rather than 
completely independent. The above result implies that 
the consensus in a real-life bird flock can also be achieved 
by the neighbor compromise mechanism, if we do not 
consider they might "copy" the flight direction of their 
common leader. 
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FIG. 3: A pair of birds (B and G) in the eleventh free flight 
show a very strong synchronization {nj — 0) and correlation 
( Cij (r) = 0.91) at the time scale from 10-270s. (a) Velocity 
envelopes (amplitudes) of B and G, and (b) phase differences 
between the pair of birds. 

The collective movement of birds, such as during the 
abrupt splitting of a flock, can instantaneously change. 
Therefore, we flrst divide a long time series into many 
short overlapping segments, and then study the pairwise 
correlation on each segment. Changing the scale of the 
time window, we can get the dynamical for multiple 
timescales. The variation of the relationship of Individ- 
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uals B and G on multiple timescales in a flight has been 
shown in Fig. 4. Because the pigeons make a circle-like 
flight in the free flights and their one-dimensional flight 
trajectory is like a pseudo-periodic time series [23, 24], 
we select the approximate cycle (20 s) as the minimum 
time window. The relationship between B and G varies 
fast with time at the small timescale [Fig. 4(a)]. Our re- 
sult shows that the collective behavior of the birds varies 
with a short timescale, and the relationship between a 
pair of pigeons is time dependent. However, With the 
timescale increasing, the frequency with which the rela- 
tionship varies reduces, for the result is an average of a 
longer time [Figs. 4(b)-(d)]. Obviously, the longer the 
time averaging is, the stabler the result of the correlation 
function becomes. 
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FIG. 4: Multi-scale pairwise correlation analysis for the differ- 
ent time- window series of B and G in the eleventh free flight. 
When the directional correlation of each time window for a 
pair of pigeons obtains its maximum value at the time delay 
Tw, T-w is called the directional correlation delay time of each 
time window. For the whole trajectories of B and G, the max- 
imum correlation coefficient C — 0.91 and the optimal time 
delay r = 0. The sampling time interval and the resolution 
of time delay both are 0.2 s. Because the pigeons make a 
circle-like flight in free flights and their one-dimensional flight 
trajectory is like a pseudo-periodic time series [23, 24], we se- 
lect the approximate cycle {w — 20 s) as the minimum time 
window. Pairwise correlation analysis for (a) the time win- 
dow is 20 s. We only calculate the optimal time delay for the 
time scale from 20-240 s (the strong synchronization segment). 
The X ax;is is the start time of the sampling window. When 
the start time of the sampling window is 20 s, Tw{t — 20.0) 
is obtained by calculating the directional correlation during 
the period from 20.0 s to 40.0 s. Next value T^{t = 20.2) is 
obtained from the period between 20.2-40.2 s, and so on up 
to the interval 220.0-240.0 s. (b) the time window is 40 s, (c) 
the time window is 100 s, and (d) the time window is 200 s. 



Recent work has illustrated the group decision rule for 
homing pigeons: compromise for small conflict and lead- 
ership for large [7] . In this study, it is difficult to measure 
the conflict level among pigeons in each flight. Neverthe- 
less, we flnd the compromise and leadership both emerge 
in the spatial domain from each dataset by calculating 
pairwise correlation [Fig. 5 and Table III]. However, we 
do not find the existence of a stable leader-follower re- 
lationship between a pair of pigeons throughout all the 
flights, and this result is different to the flndings in [7]. 
For example, two pigeons of A and B do not show the 
same relationships (different sign of t{s)) in all the flights 
[Table. II]. A leads B for three flights, A follows B for flve 
flights, and there is no leadership between them for three 
flights. That means, considering a very long time scale 
(e.g., all the flights in this study), there is no single indi- 
vidual always leading a pair of individuals. Actually, A 
leading B in this flight and then B leading A in next time 
maybe can be regarded as a general concept of reciprocal 
behaviors ( "You scratch my back, and I'll scratch yours" 
[18]). The relationship switching among group members 
supports more flexibility for individuals to response to 
external predators [25]. 



C. Coexistence of compromise and leadership 

We have shown that mutual links are ubiquitous and 
important in collective flights, so it is very valuable to in- 
vestigate whether the conclusion that certain leaders are 
able to exert more influence on the group's movement 
decisions still holds [3, 6]. We construct a flight network 
by including only those links whose maximum correla- 
tion values Cij(r) are above a given variable minimum, 
Cmin = 0.5 [Fig. 5(a)]. Here E is the follower of all the 
other individuals. Selecting different thresholds in a suit- 
able range {Cmin S [0.5 — 0.95]) to maintain the network 
connectivity gives similar results [Fig. 5(b)]. In such a 
network the nodes represent individual birds, while the 
links (arcs) denote inferred relations between their move- 
ments. If Cij(r) > Cmin and Tij > 0, we build a directed 
link pointing from the leader i to the follower j. While 
if Cij(T) > Cmin and Tij = 0, we generate a mutual link 
for the pair of pigeons. Hence we have two types of links 
in the network induced from the collective movement of 
a pigeon fiock: directed and mutual links. 

If we consider only directed links, a well-defined hier- 
archical structure is evident. On the other hand, con- 
sidering only mutual links, we find that most individuals 
share an equal structure and the individuals with recip- 
rocal links form many loops. Finally, when we take all 
the links into account, the topological structure is very 
complex because not only many directed links compose a 
hierarchical structure but also many reciprocal links form 
a large range of equal relationships. Instead of showing 
birds exclusively adopting either an averaging or a lead- 
ership strategy, our experimental analysis demonstrates 
that there is an integrated mechanism between compro- 
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TABLE II: The correlation delay time r(s) and the corresponding correlation coefficient C(r) between A and B in all the 
flights. FF means free flight and HF means homing flight. Here "— " represents that the individual does not attend the flight, 
so we have no data for analyzing its behavior. 
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0.88 
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0.95 


0.73 


0.84 


0.84 


0.99 


1.00 


0.98 



mise and leadership. 

Our data analysis shows that mutual (reciprocal) links, 
representing a coordinated correlation between a pair of 
individuals, are dense and important in pigeon flocks. A 
traditional way of quantifying the reciprocity is to com- 
pute the ratio of the number of links pointing in both di- 
rections to the total number of links i[15]: I = 
In general, the reciprocity of real-life directed networks 
ranges between the two extremes of a purely directed one 
(/ = 0, such as citation networks, where recent papers can 
cite less recent ones while the opposite cannot occur) and 
of a purely bidirectional one = 1, such as the Internet, 
where information always travels both ways along com- 
puter cables) [16]. The value of I for a real network lies 
between the above two extremes. As is shown in Ta- 
ble III, there are many mutual links and the reciprocal 
coefficient I is very high in each network of the homing 
and free flights. Moreover, the high frequency of loops 
(the size of 3) m shows that not only reciprocal relation- 
ships exist between two individuals but also such type of 
transitivity [27, 28] can be extended to three or a larger 
number of individuals. 



The number of individuals with reciprocal links and 
the entire number of individuals in a flight have been 
listed in Table III, and we find that there are very few 
pigeons having no reciprocal links. Furthermore, we list 
the ranks of individuals without reciprocal links, and we 
find such pigeons tend to be a pure leader or follower 
when observing their ranks in the group, like E in Fig. 5. 
Therefore, a reasonable conclusion is that the individ- 
ual whose position is the head or tail of the flock has a 
stronger tendency to have a directed relationship with 
others. In contrast, the individuals in the middle of the 
group tend to have a reciprocal correlation with their 
neighbors. We emphasize that the hierarchical structure 
of pigeon flocks in the previous study [6] does not im- 
ply that there arc no mutual correlations: they simply 
concentrate on the leadership aspect. Our findings show 
that both the many wrongs and leadership mechanism 
coexist in collective motion. Within and between dis- 
tinct strata, the hierarchical structure likely dominates 
long term decisions such as navigation objectives of the 
whole fiock, while the mutual interactions characterize 
local behavioral rules that are essential to maintain flock 
cohesion and alignment. 




FIG. 5: The network topology of the second free flight. If the 
directional correlation time delay Tij > and the maximum 
correlation coefficient Cij{T) > Cmin, we establish a directed 
link for the network. The gray arcs represent directed links 
from leaders to followers. If Tij = and Cij{T) > Cmin, we 
build a mutual link. The blue lines are mutual links, which 
represent a coordinated correlation between a pair of equal 
pigeons, (a) We construct the network by including only those 
edges whose maximum correlation values Cij(r) are above 
a given variable minimum, Cmin = 0.5. E is the follower 
of all the other individuals, (b) Selecting other thresholds 
(such as Cmin = 0.9) to maintain the network connectivity 
throughout the pigeon flock, also gives similar results, (c) 
The disconnected network is built by choosing Cmin = 0.99. 
Maintaining the network connectivity throughout the flight 
of the pigeon group is a key issue [26]. Therefore, it is not 
suitable to select a very high threshold to let the network lose 
the connectivity for the pigeon group. 



III. CONCLUSION AND DISCUSSION 

Our study indicates the balance between compromise 
and leadership for the organized flight of pigeons. Recip- 
rocal links represent a mutual correlation between a pair 
of individuals. Note that our result is not contradictory 
to the previous conclusion that there is a hierarchical 
structure in pigeon flocks [6], instead it is a meaningful 
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TABLE III: Summary of the statistics in each free and lioming flight performed by subjects. is the number of mutual links, 
L is the total number of links, / is the ratio of the number of links pointing in both directions L'*^ to the total number of links 
L: I = [15], m is the number of loop of size 3, A*'^ is the number of pigeons with the reciprocal relationship. A'' is the total 
number of pigeons in the flight, D is the serial number of the individuals without reciprocal links to others, and R is the ranks 
of D in each flight. " represents that the individual does not attend the flight, so we have no data for analyzing its behavior. 
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supplement. Our results show that the many wrongs and 
leadership mechanism can coexist in a collective motion. 
Hence, the dichotomy between these two mechanisms is 
false, at least for the flocking flight of pigeons. 

Our work also has significant meaning for modeling 
the collective motion of animals. If a pigeon flock only 
has a hierarchical structure, it means that the local in- 
teraction mechanism of previous models [29, 30] may not 
be adequate for simulating the group flights of homing 
pigeons (lacking the leadership), despite the interaction 
rule being dependent on the metric distance [30-32] or 
the topological distance [33]. However, our work sug- 
gests that the local interaction relationships (including 
directed and mutual links) arc sufficient to characterize 
cohesive motion of pigeon fiocks. Our results are helpful 
to provide a comprehensive picture of collective dynamic 
behavior in animal group movement and unify both the 
interaction mechanisms observed in experimental data 
and theoretical models of coherent behavior. 

We hypothesize that the integrated mechanism be- 
tween compromise and leadership also brings more ad- 
vantage to both the individuals and the whole system, to 
cope with external perturbations (e.g., predatory threat 
and food source) . Individuals of the same species come 
together to form a group because a compact fiock has 
more advantages to react to environmental perturbations 
than separate individuals [14, 34, 35]. Any external per- 
turbation for flocking movement is likely to directly cause 
a change of velocity (direction, magnitude, or both) for 
a small subset of birds that first detect the perturbation. 
Such localized changes can be transmitted to the whole 
flock to produce a collective response as if being of one 
mind [36] , making the whole group both very flexible and 
responsive. Conversely, if there were only a hierarchical 



structure in pigeon flocks, the followers cannot transmit 
any external signals to the leaders, so the whole group 
is not so sensitive to environmental perturbations. The 
sensitivity of the whole group can evolve and strengthen 
if the group members have a local interaction mechanism 
among them [37]. In particular, reciprocal links allow in- 
dividuals to interact with their neighbors and supply a 
useful way for the followers to convey information to the 
leaders. 

It needs to be noted that calculating pairwise corre- 
lation and even the methods in [6, 7] can be regarded 
as dividing a large-scale complex system (such as a fish 
school [35, 38]) into multiple local sub-systems (a pair 
of subjects). It is difficult to accurately determine how 
one individual is simultaneously affected by others (more 
than one individual), because the whole is greater than 
the sum of its parts [39]. An extension of this method 
based on the holism of complex systems [40] for analyzing 
these trajectory data should be developed in future. 
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